Pi 7i 


U.  S.  A  R  M  Y  II 

TRANSPORTATION  RESEARCH  COMMAND 
,  FORT  EUSTIS,  VIRGINIA  III 


T k ECO M  TECHNICAL  REPORT  63-7  1 

STUDY  OF 

DESIGN  FACTORS  IN  AIR  DELIVERY  FOR 
CV-7  CARIBOU  AIRCRAFT 


Task  1D643324D59806 
(Formerly  Task  9R87- 14-007-06) 

December  1963 


jWCTWBT'a—— — g-li1  UHJH  ■  -* 

-  -v 


4HJ^i'5riMU:iiJU‘-»-,'W»JP'|I"*B1*1 


iJiiJ.4li(!JLMNa| 


DISCLAIMER  NOTICE 


When  Government  drawings,  specifications,  or  other  data  are 

used  for  any  purpose  other  than  in  connection  with  a  definitely 

related  Government  procurement  operation,  the  United  States 

Government  thereby  incurs  no  responsibility  nor  any  obligation 

whatsoever;  and  the  fact  that  the  Government  may  have  formu-  *  5 

lated,  furnished,  or  in  any  way  supplied  the  said  drawings, 

specifications,  or  other  data  is  not  to  be  regarded  by  implication 

or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other  fr 

person  or  corporation,  or  conveying  any  rights  or  permission, 

to  manufacture,  use,  or  sell  any  patented  invention  that  may  in 

any  way  be  related  thereto. 


*  *  * 

DDC  AVAILABILITY  NOTICE 

Qualified  requesters  may  obtain  copies  of  this  report  from 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 


*  *  * 

This  report  has  been  released  to  the  Office  of  Technical  Services, 
U.  S.  Department  of  Commerce,  Washington  25,  D.  C.  ,  for  sale 
to  the  general  public. 


*  ★  $ 

The  findings  and  recommendations  contained  in  this  report  are  y 

those  of  the  preparing  agency  and  do  not  necessarily  reflect  the 
views  of  the  U.  S.  Army  Mobility  Command,  the  U.  S.  Army 
Materiel  Command,  or  the  Department  of  the  Army. 


Task  1D643324D59806 
(Formerly  Task  9R87-14-007-06) 

TRECOM  Technical  Report  63-71 

December  1963 

I 

STUDY  or 

DESIGN  FACTORS  IN  AIR  DELIVERY  FOR 
CV-7  CARIBOU  AIRCRAFT 


Prepared  by 
R.  D.  H.  Jones 


APPROVED. 


U.  S.  ARMY  TRANSPORTATION  RESEARCH  COMMAND 
FORT  EUSTIS,  VIRGINIA 


CONTENTS 


Page 

LIST  OF  SYMBOLS .  iv 

I 

SUMMARY .  1 

CONCLUSIONS .  3 

INTRODUCTION .  4 

ASSUMPTIONS .  5 

ONE-DEGREE-OF-FREEDOM  MOTION .  7 

TWO-DEGREES-OF-FREEDOM  MOTION .  11 

THREE-DEGREES-OF-FREEDOM  MOTION  ....  18 

APPENDIX.  FORTRAN  PROGRAMS .  25 


DISTRIBUTION  . 


29 


SYMBOLS 


Alphabetical 

Character 

AM 

P 

DE  or  DD 

8 

C 

H 

HG 

HH 

*G 

KE 

L 


fortran 

Character 

AM 

S 

DE  or  DD 


HG 


HH 


Q 


QD 


length  of  container 
Tan' 


r-lf2HG\ 
V  AM J 


arbitrary  interval  of  time  between 
successive  points  in  solution  of  a 
differential  equation 

acceleration  of  center  of  gravity 

center  of  gravity 

mid -point  of  platform  of  container 

height  oi  center  of  gravity  of  con¬ 
tainer  above  underside 

height  of  container 

moment  of  inertia  of  container  about 
center  of  gravity,  about  horizontal 
axis  perpendicular  to  line  of  flight 

kinetic  energy 

sill  ol  floor 

force  oi  extraction  of  parachute  on 
container  in  gravitational  units 

distance  of  mid-point  of  underside 
of  container  from  aft  edge  of  floor 

first  derivative  of  Q  with  respect 
to  time 


iv 


Alphabetical 

Character 

FORTRAN 

Character 

t 

time 

V 

potential  energy 

W 

weight  of  container  in  gravitational 
units 

f 

* 

X 

X 

distance  in  inches  from  forward 
bulkhead  of  cargo  compartment  to 
some  indicated  point  of  container 

xst 

XST 

distance  from  bulkhead  of  point  on 
floor  from  which  container  starts 
from  rest 

XB^B 

coordinates  of  corner  B  of  container 

Vc 

coordinates  of  center  of  gravity  of 
container 

e 

T  or  TH 

angle  through  which  container  has 
turned  from  original  horizontal 
position 

•  • 

e 

TD 

first  derivative  of  B  with  respect  to 
time 

0 

TDD 

second  derivative  of  0  with  respect 
to  time 

X 

YAL 

angle  of  inclination  of  parachute 
cord  to  horizontal 

* 

<P 

F 

0  +  X  +  8 

* 

• 

FD 

first  derivative  of  <f>  with  respect 
to  time 

FDD 

second  derivative  of  <J>  with  respect 
to  time 
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SUMMARY 


A  study  is  made  of  the  motion  of  containers  loaded  in  the  cargo  com¬ 
partment  of  the  Caribou  CV-7  aircraft  when  the  containers  are  to  be 
extracted  by  parachute.  *  The  cord  of  the  parachute  is  attached  to  the 
foot  of  the  aft  side  of  the  container  (a  method  that  gives  the  container  a 
distinct  path  of  descent)  rather  than  to  the  mid-point  of  the  aft  side. 
None  of  the  containers  loaded  in  any  of  the  positions  defined  in  the  six 
standard  loading  diagrams  used  in  this  report  will  strike  either  the 
ceiling  or  the  ramp  at  the  rear  of  the  cargo  compartment.  The  safe 
envelope  for  any  load  is  a  rectangle  of  76  inches  in  height  above  the 
floor;  this  will  have  to  be  further  reduced  to  allow  for  the  height  of  the 
rollers  above  the  floor. 

The  motion  of  the  container  from  the  instant  when  it  is  released  from  a 
stationary  position  on  the  floor  of  the  aircraft  and  is  allowed  to  move 
aft  under  the  force  of  the  extraction  parachute  until  it  is  falling  freely 
under  its  main  chute  can  be  divided  into  the  following  three  phases:** 

One-Degree-of-Freedom  Motion 

During  this  phase,  the  container  slides  down  the  rollers.  This 
sliding,  with  the  underside  of  the  container  resting  horizontally  on 
the  rollers,  continues  for  a  short  distance  after  the  mid-point  of 
AB  on  the  container  has  reached  point  L.  While  the  108-,  162-, 
and  216-inc’n-long  containers  slide  normally  down  the  floor,  the 
54-inch-long  container  (see  Figure  1)  possesses  a  special  motion 
of  its  own.  Provided  it  is  clear  of  the  guide  rail,  comer  A  rises 
very  slightly  above  the  floor  but  returns  at  once  to  the  horizontal 
position. 


*  A  subsequent  report  will  discuss  the  minimum  height  from  which  containers  may  be  dropped  in  order  to  strike 
the  ground  when  the  parachute  is  located  directly  above  the  container 

**  A  rigid  body  moving  in  space  possesses  six  degrees  of  freedom,  but  we  consider  that  the  fore  and  aft  vertical 
sides  of  the  container  remain  throughout  in  the  same  vertical  planes.  A  fixed  plane  *r  the  co.  ta<ner  is  therefore 
constrained  to  pass  through  three  arbitrary  points,  and  in  this  way  three  of  the  six  degrees  rf  fxaedom  are 
"ankylosed".  (Reference  J.  H.,  Jeans,  The  Dynamical  Theory  of  Gases,  4th  Edition,  Cambridge  Jmveraity  Press f/ 
New  York,  New  York,  19Z5.  p.  67,,  par.  77.) 


Direction  of  Motion  of  Aircraft 


Two -Degrees  -of -Freedom 
Motion 

The  container  continues  to 
move  aft  with  its  underside 
in  contact  with  the  aft  edge 
of  the  floor,  L.  The  initial 
conditions  for  this  two -degree 
motion  are  the  final  conditions 
for  the  one-degree  motion. 

The  container  always  remains 
in  contact  with  edge  L.  The 
possibility  that  corner  C  of 
the  container  might  come 

dangerously  close  to  the  ceiling  of  the  cargo  compartment  must 
be  watched.  The  ceiling  is  82  inches  above  the  floor;  if  a  clear¬ 
ance  of  6  inches  is  allowed,  a  safe  envelope  to  study  for  a  container 
or  other  load  is  one  that  can  have  no  point  that  would  foul  a  height 
76  inches  above  floor  level;  here,  again,  it  is  necessary  to  make 
allowance  for  height  of  rollers. 


Figure  1. 


Movement  of  Container  in 
Aircraft  With  One-Degree- 
of-Freedom  Motion. 


Three-Degrees-of-Freedom  Motion 

This  phase  occurs  after  B  has  cleared  L  and  the  container  is  falling 
freely  under  the  force  of  gravity  and  the  force  of  its  main  chute. 

The  initial  conditions  for  this  three-degree  motion  are  provided  by 
the  final  conditions  for  the  two-degree  motion. 

When  data  were  computed  for  this  report,  it  was  assumed  that  the  ramp 
was  inclined  to  the  horizontal  at  an  angle  of  15  degrees,  mat  there  was 
no  friction  between  the  containers  and  the  floor,  and  that  the  force 
exerted  by  the  extraction  parachute  is  1.  5  times  that  of  gravity  at  all 
stages  of  the  process  of  extraction. 

The  analytical  mechanics  approach  to  the  trajectories  of  the  containers 
is  described.  The  work  was  performed  on  an  IBM  1620  computer,  and 
the  more  important  FORTRAN  programs  used  to  obtain  the  data  are 
presented  in  the  appendix. 
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CONCLUSIONS 


mmm** 


It  is  concluded  that: 

1.  None  of  the  containers  loaded  in  any  of  the  positions  shown  in 
the  six  loading  plans  will  hit  the  ceiling  or  the  ramp  when  the 
latter  is  inclined  at  15  degrees  to  the  horizontal. 

i 

2.  The  safe  envelope  parallel  to  the  axis  of  the  aircraft  is  a 
rectangle  of  75  inches  in  height  above  the  floor. 

3.  Small  variations  in  the  angle  of  inclination  of  the  parachute 
cord  to  the  horizontal  do  not  affect  these  results. 

4.  Since  the  containers  slide  horizontally  down  the  floor,  the 
length  of  the  guide  rail  has  no  effect  on  the  motion. 


i 
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INTRODUCTION 


This  report  is  written  to  supplement  "A  Study  of  Air  Delivery  Supply 
Systems  for  the  CV-7  Caribou  Aircraft",  published  by  the  Quarter¬ 
master  Research  and  Engineering  Center,  Aerial  Delivery  Systems 
Office,  Natick,  Massachusetts,  in  June  I960. 


A  report  prepared  in  May  1963  by  the  U.  S.  Army  Transportation  Re¬ 
search  Command  (USATRECOM)  showed  an  approximate  method  for 
calculating  the  trajectories  of  containers  to  be  extracted  by  parachute 
from  various  stationary  positions  on  the  aircraft  floor  when  the  para¬ 
chute  cord  was  attached  to  the  mid-point  of  the  aft  side  of  the  container. 
Attaching  the  cord  to  the  mid-point  gives  a  distinctly  different  motion 
to  the  container.  {The  higher  the  point  of  attachment,  the  greater  will 
be  the  tendency  of  the  container  to  tip  and  strike  the  ceiling.)  Gener¬ 
alized  equations  of  motion  were  not  employed  in  the  report,  and  the 
trajectories  obtained  were  not  exact.  The  results  indicated  that  con¬ 
tainers  starting  from  positions  300  inches  or  more  from  the  forward 
bulkhead  would  strike  the  ramp. 

The  present  report  deals  with  the  trajectory  of  the  container  when  the 
cord  is  attached  to  the  foot  of  the  aft  side;  attachment  at  this  point 
gives  the  container  a  totally  distinct  path, 

NOTE:  Attaching  the  cord  at  an  intermediate  point  between  the 

foot  and  mid-point  of  the  aft  side  could  also  be  considered. 


Information  was  required  regarding  the  effect  of  the  length  of  the  guide 
rail  on  the  container  trajectory.  Theoretically,  a  container  extracted 
by  a  cord  attached  to  the  foot  of  the  aft  side  remains  virtually  hori¬ 
zontal  until  it  reaches  the  aft  end  of  the  floor  of  the  cargo  compartment. 
Hence,  on  the  model  considered  for  this  report,  the  length  of  the 
guide  rail  has  little  or  no  effect  on  the  container  path.  However,  to 
load  containers  into  the  aircraft,  the  top  of  the  guide  rail  must  stop 
short  of  the  end  ot  the  platform  by  a  distance  approximately  equal  to 
half  the  length  of  the  longest  container.  The  longest  container  is  216 
inches,  so  the  horizontal  flange  of  the  rail  must  stop  around  108  inches 
short  of  the  rear  end  of  the  floor. 
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Since  the  container  remains  horizontal  while  sliding  down  the  floor 
whether  there  is  a  guide  rail  or  not,  the  possibility  of  the  roof's  being 
struck  by  the  container  during  this  first  phase  of  the  motion  need  not 
be  envisaged.  A  FORTRAN  program  has  been  developed  to  find  the 
maximum  safe  envelope,  but  the  absence  of  appreciable  tipping  means 
that  a  container  76  inches  in  height  above  the  floor  will  not  strike  the 
ceiling;  allowance  must  be  made  for  the  height  of  the  rollers. 

Six  loading  plans  for  containers  of  four  different  lengths  in  the  cargo 
compartment  were  considered  in  this  study  and  are  shown  in  Figure  2; 
the  shaded  areas  represent  containers  of  the  lengths  indicated,  loaded 
horizontally. 


ASSUMPTIONS 


It  has  been  assumed  that:. 

1.  If  the  tail  ramp  is  inclined  at  15  degrees  to  the  horizontal,  no 
container  considered  in  the  six  loading  plans  will  strike  the 
ramp. 

2.  Unless  otherwise  stated,  the  cord  of  the  extraction  parachute 
remains  horizontal.  Two  of  the  four  FORTRAN  programs 
presented  in  the  appendix  include  as  a  parameter  the  angle  of 
inclination  of  the  cord,  and  in  these  programs  a  non-zero 
value  can  be  given  to  this  angle  without  difficulty.  * 

3.  The  force  exerted  by  the  extraction  parachute  is  1.5  times 
that  of  gravity  at  all  stages  of  the  process  of  extraction.  This 
is  the  sor+  of  practical  factor  regarding  which  experimental 
data  would  be  valuable. 

4.  There  is  no  friction  between  the  container  and  the  floor.  (It 
is  natural  to  wonder  what  effect  the  presence  of  some  friction 
would  make.  Friction  would  reduce  the  effective  force  of  the 
chute  cord  below  1. 5  times  that  of  gravity;  by  so  doing,  it 
would  decrease  the  horizontal  velocity  acquired  by  the  con¬ 
tainer  while  sliding  down  the  floor  and  thus  increase  the 
chance  of  the  ramp's  being  struck.  However,  if  the  cord 


*  Since  four  principal  FORTRAN  programs  have  been  included  m  this  report,  agencies  with  computing  facilities 
will  be  able  to  extend  the  calculations  to  further  ranges  of  the  parameters.  However,  USATRECOM  will  be  glad 
to  amplify  any  of  the  calculations  for  those  who  cannot  obtain  the  use  of  a  computer. 
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Loading  Plans  for  Cargo  Compartment  of  Caribou  CV-7  Aircraft. 


remained  attached  to  the  foot  of  the  aft  side  of  the  container, 
there  would  be  no  tendency  to  tip,  and  thus  the  possibility  of 
the  ceilings's  being  hit  would  not  be  increased. ) 


ONE-DEGREE-OF-FREEDOM  MOTION 


t 


t 


1 


i 


1 
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Although  the  first  phase  of  the  motion  has  been  referred  to  as  having 
one  degree  of  freedom,  the  configuration  of  the  container  in  Figure  1 
is  described  in  terms  of  the  two  coordinates  x  and  0.  It  will  be  shown 
that  0  departs  from  zero  value  for  a  short  time  only.  In  effect,  there¬ 
fore,  this  phase  is  one  of  a  single  degree  of  freedom.  Let  x^,  and  y^ 

be  the  coordinates  of  G,  the  center  of  gravity,  with  respect  to  the  axes 
shown;  note  that  the  y  axis  has  positive  direction  upward  in  this  case. 


If 

and 

then 

and 


xG  *  x  +  y  cos  9  -  HG  sin  0 

AM 

yG  *  Y  8*n  9  +  HG  cos  0, 

•  •  •  All 

Xg  *  x  -  0<Y  sin  0  +  HG  cos  0) 

yr  =  0(4r  cos  0  -  HG  sin  0). 

G  2 


The  kinetic  energy  of  the  container  is  equal  to  that  of  its  total  mass 
moving  with  its  center  of  gravity  plus  the  kinetic  energy  of  rotation 
about  the  center  of  gravity.  If  it  is  assumed  that  the  container  is 
symmetrical  in  its  loading  and  construction,  then  the  height  equals 
twice  HG,  and  the  radius  of  gyration  about  a  horizontal  axis  through  G 
perpendicular  to  the  plane  of  the  paper  is 


AM2  +  4HG2 
12 


If  KE  is  the  kinetic  energy,  then 


2KSS  .  »Ci»  ♦  .  ^-VHGS]  . 

Substituting  for  x  and  y,  this  becomes 

2KEg  «  W[x2  -  2x0(&5  sin  0  +  HG  coa  0)  ♦  02(^— . 

2  ° 
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As  there  is  no  friction  between  the  container  and  the  floor,  the  con¬ 
tainer  has  a  potential  energy  function  under  the  action  of  the  forces  of 
gravity  and  the  pull  of  the  chute.  This  potential  function  is 

AM 

V  s  constant  +  #(-3-  sin  0  +  HG  cos  6)  -  P(x  ♦  AM  cos  0),. 

LAGRANGE  EQUATION  OF  MOTION  FOR  GENERALIZED  COORDI¬ 
NATES  x  AND  6 

The  Lagrange  equation  of  motion  for  6  is 

i sin  0  -  xHG  cos  0  +  9( AM2-y^-G-l)]  -  cos  0  -  HG  sin  6)] 

5  -Wfffl  cos  0  -  HG  sin  0)  -  P  AM  sin  0  . 

When  the  differentiation  is  performed,  this  equation  is  reduced  to 

-x(—  sin  0  +■  HG  cos  0)  -  x0(—  cos  0  -  HG  sin  0)  +  9  (AM?-  f„„fflg£) 

2  2  3 

+x6 ( cos  0  -  HG  sin  0)  =  g(-  cos  0  +  HG  sin  0  -  1.5AM  sin  0) 

P 

on  dividing  through  by  W  and  setting  —  equal  to  1.5  in  the  last  term. 

w 

•  • 

The  product  term  in  x0  cancels  out,  and  the  equation  becomes 


-x(AE  sin  0  +  HG  cos  0)  +  ±  HHS2) 

2  3 


=  g( HG  sin  0  -  y  cos  0  -  1,5AM  sin  0).  (1) 


The  Lagrange  equation  of  motion  corresponding  to  the  coordinate  x  is 

t  -Ifx  -  0(-~  sin  0  +  HG  cos  0)3  =  P. 
g  df  ? 

P 

When  the  differentiation  is  performed  and  ^  is  set  equal  to  1.5,  the 
equation  becomes 

x  -  0(M  sin  e  *  HG  cos  0)  -  02(^i  cos  0  -  HG  sin  0)  »  1.5g  .  (2) 


i 


r 


8 


Whether  the  initial  value  of  8  is  positive  or  negative  must  be  deter¬ 
mined.  When  this  phase  of  the  motion  starts,  whether  the  container  is 
held  horizontal  by  the  the  guide  rail  or  not,  both  6  and  8  are  zero.  If 
6  and  o  are  set  equal  to  zero  in  Lagrange  equations  (1)  and  (2),  then 

-xHG  ♦  5(«i?  t  .  ,«() 

3  2 

and 


x  -  0HG  =  1.5g  . 


Hence, 


g 


■HG 


AM 

"  T 


-HG 


AM2  +  4HG2 
3 


1  1.5 


1  -HG 


The  denominator  of  the  right-hand  member  is  always  negative.  There¬ 
fore,  if  3HG  is  greater  than  AM,  8  is  initially  positive.  AM  =  54  inches 
for  the  small  container,  and  this  condition  is  always  satisfied.  Hence, 

6  is  initially  positive  for  the  small  container;  that  is,  it  tips  backward 
as  in  Figure  1.  However,  when  the  two  Lagrange  equations  are  solved 
numerically,  8  at  once  returns  to  zero  value.  For  the  other  three 
sizes  of  the  container,  AM  is  greater  than  3HG;  as  a  result,  no  tipping 
occurs,  and  the  container  pursues  a  true  one-degree-of-freedom 
motion  until  it  reaches  L,  the  end  of  the  platform,  and  begins  to  tip  as 
explained  in  the  section  on  two-degrees-of-freedom  motion.  In  the 
latter  case,  the  tipping  is  in  the  opposite  direction;  that  is,  corner  A 
of  the  container  is  not  above  corner  B  as  shown  in  Figure  1  but  is  below 
B. 


FORTRAN  NUMERICAL  SOLUTION  OF  EQUATIONS  (1)  AND  (2) 

Since  the  initial  value  of  8  has  been  discussed  and  the  criterion  for  its 
being  positive  has  been  observed,  the  solving  of  equations  (1)  and  (2) 
can  be  continued  in  order  to  determine  x  and  9  for  subsequent  values 
of  time.  FORTRAN  Program  Number  Two  (see  appendix)  provides  an 
approximate  numerical  solution  of  the  two  ordinary  nonlinear  differen¬ 
tial  equations  in  x  and  8. 


The  values  of  the  parameters  AM  and  HG  are  read  into  the  computer. 

On  the  same  line  of  the  program,  DE  is  also  read  in.  (This  is  an 
arbitrary  interval  of  time  which  was  taken  as  .01  second  between  suc¬ 
cessive  configurations  of  the  system  at  which  solutions  of  the  equations 
were  sought. )  In  the  next  line  of  the  program,  the  initial  value  of  x, 
called  X(l),  is  read  in;  this  is  the  starting  station  of  the  forward  side 
of  the  container.  In  the  same  line  of  the  program,  the  following  is 
read  in: 

XD(1)  the  first  derivative  of  x  with  respect  to  time  at  the  initial 
instant; 

T(l)  the  initial  value  of  6;  and 
TD(1)  the  initial  value  of 

XD(1),  T(l),  and  TD(1)  all  have  the  value  of  zero  in  the  case  of  a 
small  container  starting  from  the  rearmost  station.  If,  however,  the 
container  were  released  under  the  guide  rail,  X(l)  would  be  the  length 
of  the  guide  rail,  and  XD(1)  would  be  the  velocity  the  container  ac¬ 
quires  in  traveling  down  the  guide  rail- -that  is,  from  whatever  station 
under  the  guide  rail  it  started  from  rest. 

Lines  5  through  12  of  FORTRAN  Program  Number  Two  give  the  second 
differential  coefficients  of  x  and  of  0,  that  is,  XDD(  1)  and  TDD(l). 

Lines  14  through  17  are  short  Taylor-series  equations  to  find  the  valuer 
of  X  and  T  and  of  their  first  derivatives  at  the  close  of  the  next  time 
interval.  Lines  18  through  26  form  and  solve  two  ordinary  simultaneous 
linear  algebraic  equations  for  TDD  and  XDD,  which  in  turn  are  sub¬ 
stituted  in  the  Taylor  series  of  lines  14  through  17;  the  whole  process 
is  repeated  as  often  as  desired. 

If  AM  -  54  inches,  HG  =  30  inches,  DD  =  .01  second,  and  X(l)  =  291 
inches,  the  successive  values  for  T,  or  0,  are  as  follows: 

T(l)  =  .000000  radian  T(3)  =  .000367  radian 

T(2)  =  .000273  radian  T(4)  =  000444  radian 

A  negative  value  of  0  is,  of  course,  impossible  as  it  is  constrained  to 
be  non-negative  by  the  platform.  In  physical  terms,  this  means  that  tht 
small  container  makes  an  initial  jerk  with  positive  0  as  shown  above 
for  T(l)  to  T(3)  but  returns  to  the  horizontal  position  and  slides  down 
the  floor. 

If  the  length  of  the  container  were  reduced  so  that  AM  became,  say,  1 
inch  only,  then  the  force  P  would  tip  it  over  completely.  Somewhere 
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between  the  54-inch  and  the  1-inch  containers,  there  must  be  some 
value  of  AMq  which  divides  small  containers  into  those  that  tip  slightly 
and  return  to  the  horizontal  and  those  that  fall  completely  over.  Such 
a  value  for  AM„  would  mark  what  Loupounoff  called  a  "point  of  dynami¬ 
cal  bifurcation". 

Thus,  for  all  four  types  of  containers,  the  first  phase  of  the  motion  is, 
in  effect,  a  sliding  down  the  floor  until  the  end,  L,  is  reached.  For 
this  reason,  the  first  phase  is  referred  to  as  the  one-degree-of- 
freedom  phase. 


TWO- DEGREES -OF- FREEDOM  MOTION 


For  this  phase  of  the  motion,  0  (the  angle  of  inclination  of  the  con¬ 
tainer  to  the  horizontal  reckoned  positive  as  shown  in  Figure  3)  and  q 
(the  distance  from  L  to  the  mid-point  of  the  underside,  AB,  of  the  con¬ 
tainer)  are  selected  as  general  coordinates. 


Let  xS£  be  the  distance  in 
inches  between  the  bulkhead 
at  the  forward  end  of  the 
cargo  compartment  and  the 
forward  end  of  the  container 
when  the  latter  starts  to 
move  at  the  beginning  of  the 
one-degree  phase.  Let  Q(l) 
be  the  value  of  q  when  tipping 
starts.  Then,  between  re¬ 
lease  from  station  xs^  and 
start  of  tipping,  the  container 
has  moved 

(377  -  xst  +  Q(l)  -  1/2  length  of  container)  inches 

2 

under  an  acceleration  of  1.5g  inches  per  second  .  Hence,  the  final 
velocity  of  the  container  at  the  close  of  phase  one  and  at  the  start  of 
phase  two  is  6. 

V3g[377  -  x  TWToUri  inches  per  second 
st  2 

*"^-^9.2[377  -  xat  -  +  Q(1)T  inches  per  second  (3) 
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Direction  of  Motion  of  Aircraft 
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Figure  3.  Sketch  Illustrating  Two- 
Degrees -of-Freedom 
Motion 
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For  phase  two,  this  is  the  initial  value  of  the  first  derivative  of  q  and 
is  called  QD(1).  The  initial  values  of  6  and  6  are  both  zero;  these  are 
referred  to  as  T(l)  and  TD(1)  respectively. 

Taking  axes  of  coordinates  as  shown  in  Figure  3  (that  is,  with  positive 
direction  of  y  axis  downward),  the  center  of  gravity  of  the  container,  G, 
has  coordinates 

Xq  =  q  cos  0  +  HG  sin  0  and  »  q  sin  0  -  HG  cos  0. 

The  kinetic  energy  of  the  container  is  equal  to  that  of  its  total  mass 
moving  with  its  center  of  gravity  together  with  the  kinetic  energy  of 
rotation  about  its  center  of  gravity.  Therefore,  if  KE  is  its  total 
kinetic  energy, 

2KEg  *  W(i|  ♦  y2)  +  IQ02 

s  W([q  cos  0  ♦  0(-q  sin  0  HG  cos  0)]2 

*  +  [q  sin  0  ♦  0(q  cos  0  +  H6  sin  O)]2}  +  I„02 

fa 

•  WCq2  ♦  02(q2  +  HG2)  t  2HGq0]  ♦  1Q02. 

Since  the  motion  is  assumed  to  be  frictionless,  the  container  has  a 
potential  energy  function  under  the  force  of  gravity  and  the  force  of 
j  the  chute;  thus, 

i 

1  V  =  constant  +  g[-W(q  sin  0  -  HG  cos  0) 


-  P  cos  A(q  +  -M)  cos  0  ♦  P  sin  A(q  ♦  — )  sin  0] 

2  2 


8  constant  ♦  g[-W(q  sin  0  -  HG  cos  0)  -  P(q  +  AS)  cos  (0  +  A)]. 


THE  LAGRANGE  EQUATIONS  OF  MOTION  FOR  0  AND  q 
The  Lagrange  equation  of  motion  for  0  is 
^[IgO  +  W(q2  +  HG2)  0  +  HG  W  q] 

s  W  a(  n  f'na  ft  »  HC  etn  A)  —  Pfn  + 

2 
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=  Wg(q  cos  0  ♦  HG  sin  0)  -  P(q  t  AS.)  g  sin  (0  ♦  A)t 


- - — — ■******"&&—  iffiyyaiiMl 
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and  the  equation  for  q  is 


W-^<q  ♦  0HG)  -  W02q  »  g[W  sin  6  ♦  P  cos  (0  ♦  X)]. 


Divide  these  equations  through  by  W,  and  let  ^  =1.5.  Also,  let 


.-9  AM2  +  4HG2 
G  *  Ig  ♦  HG2  *  . jy- - ♦  HG2, 


using  the  known  formula  for  radius  of  gyration  of  a  rectangular  paral¬ 
lelepiped  about  an  axis  perpendicular  to  its  plane  and  through  its  cen¬ 
ter  of  gravity.  Thus, 

6(G  +  q2)  +  2qq0  +  HGq 

AH 

*  g[q  cos  0  ♦  HG  sin  0  -  1.5(q  ♦  -5-)  sin  (0  +  X)]  (4) 

*  ’  - 

and 

HG0  -  02q  ♦  q  «  g[sin  0  ♦  1,5  cos  (0  ♦  X)]  . 


FORTRAN  NUMERICAL  SOLUTION  OF  EQUATIONS  (4)  AND  (5) 

The  method  for  solving  these  equations  by  using  FORTRAN  Program 
Number  One  is  as  follows: 

d6  dq 

First,  calculate  the  initial  values  of  0,  q,  ^7-,  and  in  the  program 

these  symbols  are  called  T(l),  Q(l),  TD(1),  and  QD(1).  Substitute 

these  values  in  the  two  differential  equations,  and  solve  them  for  0  and 
•• 

q  as  two  ordinary  simultaneous  linear  algebraic  equations  in  these  two 
unknowns . 

d0 

Second,  by  means  of  Taylor's  theorem,  find  the  values  of  0,  q,  jp  and 
at  .01  second  after  the  initial  instant.  Thus, 

01  2 

T(  2)  =  T(  1)  +  .01  TD(  1)  +  (-^i)  TDD(l) 


and 


TD(2)  =  TD(1)  ♦  .01  TDD(l)  . 

Similarly,  find  the  values  for  Q(2)  and  QD(2). 

Third,  substitute  these  values  of  0,  q,  and  ^  in  the  two  differen¬ 
tial  equations  and  solve  them  as  simultaneous  algebraic  equations  for 
tne  second  derivatives  of  0  and  q. 


Repeat  the  same  process  for  T(3)  and  Q(3),  etc. 


Calculation  of  QD(1) 

dq 

Instead  of  reading  in  the  value  of  “ft"  at  the  initial  instant  as  was  done 

before,  in  certain  of  the  FORTRAN  programs  the  following  line,  which 
calculates  QD(1)  from  equation  (3),  is  inserted  in  the  program: 


QD(1)  s  V1159.2C377  -  x^ 


+  Q(l)]. 


Initial  Value  of  q 


The  calculation  of  the  initial  value  q[that  is,  Q{1)  ]  possesses  some 
interesting  features.  If  it  is  assumed  that  this  initial  value  is  zero 
(that  is,  that  tipping  starts  when  the  CciUCii  Ol  gravity  of  the  container 
is  vertically  ovor  point  L.  a l  the  end  of  the  floor),  then  the  initial  value 

4-'0 

of  ' — -  is  negative,  and  hence,  0  starts  to  assume  negative  values.  In 
dt  1 

physical  terms,  thi3  means  that  the  container  begins  to  tip  into  the 
floor.  Thus,  coordinates  0  and  q  must  be  regarded  as  being  subject  to 
a  constraint:  namely,  0  must  be  positive.  If  0  and  are  equal  to 

zsrc  in  the  differential  equations, 

0(G  +  q2)  +  HGq  =  g[q  -  1.5(q  t  4jr)  sin  X] 


q  =  1.5g  cos  X. 


Solving  for  9, 


“  q(  1  -  1.5  sin  X)  -  1.5(4p  sin  X  -  1.5HG  cos  X 
g  G  +  q2  -  HG2  ~ 


The  denominator  of  the  above  is  always  positive;  hence,  for  9'  to  be 
positive  initially, 

q(l  -  1.5  sin  X)  must  be  >  1.5(HG  cos  X  ♦  sin  X). 

The  sign  of  equality  in  the  above  equation  gives  the  initial  value  of  q 
and  is  incorporated  in  certain  of  the  programs  to  avoid  readme  in  Q(l) 
as  data.  It  is  evident  that  in  no  case  can  Q(l)  be  greater  than -j—  . 


Reaction  of  Aft  End  of  Platform  on  Container 


The  fundamental  equations  for  the  two-degrees-of-freedom  motion 
have  been  developed  on  the  assumption  that  from  the  moment  that 
corner  A  (see  Figure  3)  of  the  container  reaches  L,  the  container  re¬ 
mains  in  contact  with  the  aft  edge  of  the  cargo  compartment  until 
corner  B  of  the  container  reaches  L.  In  order  to  make  this  assumption, 
the  force  of  reaction  of  the  edge  L  on  the  container  must  be  considered 
positive.  Since  the  container  has  a  smooth  lower  surface,  the  direc¬ 
tion  along  which  the  reaction  acts  is  normal  to  the  underside  of  the 
container.  Let  R  units  of  force  be  the  magnitude  of  the  reaction;  then, 
if  the  forces  are  resolved  vertically, 


(j2y 

W  -  P  sin  X  -  R  ccs  6  =  —  — —  j? 

8  dt2 

where  Yq  is  reckoned  positive  downward  as  in  Figure  3.  When  the  two 
differentiations  are  carried  out, 

R  cos  e  =  W  -  P  sin  X  +-^(-HG  cos  002  -  HG  sin  00 

o 

-  q  sin  0  -  2q  cos  00  +  q02  sin  0  -  q0  cos  0). 

By  means  of  this  relation,  the  sign  of  R  may  be  tested  for  being  posi¬ 
tive  or  negative  by  setting  up  a  short  FORTRAN  subprogram  which  can 
be  incorporated  in  the  main  dimensioned  program  which  solves  the 
two-degrees-of-freedom  Lagrange  equations  and  gives  q  and  0  and 
their  first  two  derivatives.  These  six  quantities  [Q(l),  QD(1),  QDD(l), 
T(l),  TD(1),  and  TDD{1)  }  remain  in  core  storage,  since  they  are 
dimensioned  and  are  thus  available  for  the  subprogram.  The  sub¬ 
program  shows  that  R  is  positive  in  all  cases. 

Accuracy  of  Solution  of  Differential  Equations 

A  subprogram  has  been  incorporated  in  FORTRAN  Program  Number 
One  to  test  the  accuracy  of  the  solutions.  This  subprogram  transfers 
all  the  terms  to  the  left-hand  members  of  the  equations  and  substitutes 
in  them  the  values  found  for  q  and  0  and  their  first  two  derivatives. 

This  gives  the  error;  expressed  as  a  proportion  of  the  sum  of  the 
absolute  values  of  the  terms  in  the  differential  equations,  the  errors 
are  around  one  part  in  1,000. 
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INCLINATION  OF  CONTAINER 


Table  1  presents  some  of  the  numerical  results  of  FORTRAN  Program 
Number  One  (see  appendix)  which  solves  the  equations  of  motion  for 
the  two-degree  phase.  The  line  of  data  following  the  "END"  of  this 
program  relates  to  a  216-inch-long  container  placed  with  its  center  of 
gravity  at  the  center  of  gravity  of  the  aircraft;  that  is,  with  the  forward 
side  of  the  container  77.  8  inches  from  the  forward  bulkhead.  The 
height  of  the  center  of  gravity  of  the  container  in  this  case  is  considered 
to  be  20  inches — a  reasonable  value  for  a  load  such  as  a  small  truck. 

TABLE  I 

TWO-DEGREES-OF-FREEDOM  MOTION  FOR  216-INCH-LONG  CONTAINER 


Time  From  Start 

Angle  of  Inclination 

Value  of  q 

of  Tipping 

to  Horizontal 

(seconds) 

(radi?  ns) 

(inches) 

.00 

.000000000 

37.7 

.01 

.000000323 

42.9 

.02 

.0000174 

48.1 

.03 

. 0000800 

53.4 

.04 

.000210 

58.8 

.05 

.000423 

64.2 

.06 

. 000728 

69.6 

.07 

.00113 

75.2 

.08 

.00164 

80.7 

.09 

. 00225 

86.4 

.  10 

.00295 

92.0 

.  11 

.00375 

97.8 

.  12 

. 00464 

103.6 

.  13 

. 00562 

109.4 

No  further  entries  have  been  made  in  the  table  because  q  is  already 
greater  than  half  the  length  of  the  container;  consequently,  the  container 
has  cleared  the  end  L  of  the  floor.  The  angle  of  inclination  of  the  para¬ 
chute  cord  to  the  horizontal  has  been  taken  as  .  034907  radian,  that  is, 
two  degrees  of  arc.  The  table  shows  clearly  that  the  inclination  of  the 
container  is  slight  under  representative  conditions. 


CALCULATION  OF  SAFE  ENVELOPE 


It  has  been  shown  that  a  container  tips  through  so  small  an  angle  that  it 
seems  almost  unnecessary  to  develop  a  FORTRAN  program  to  find  the 
safe  envelope.  It  could  be  assumed,  for  example,  that  a  76-inch-high 
container  would  at  no  time  during  the  two-degree  phase  have  any  point 
more  than  77  inches  above  the  floor  and  hence  that  a  constant  height  of 
76  inches  would  be  safe  in  view  of  the  6-inch  clearance  provided.  This 
is  correct  for  all  the  containers  considered  in  the  six  loading  plans. 

FORTRAN  PROGRAM  NUMBER  FOUR 

In  Figure  4,  let  AB  be  the 
underside  of  the  container, 

H  the  mid-point  of  this 
underside,  and  q  the  dis¬ 
tance  from  H  to  L.  Let  NU 
be  an  ordinate  of  the  con¬ 
tainer  (that  is,  NU  is  per¬ 
pendicular  to  AB).  In  lines 
1  through  34  of  FORTRAN 
Program  Number  Four  (see 
appendix),  the  values  of  q 
and  0  are  found  as  previously.  Figure  4.  Sketch  for  Determining 
In  lines  36  and  37,  N  is  set  Safe  Envelope, 

at  the  end  A  and  successively 

at  4-inch- intervals  along  the  underside,  AB,  of  the  container.  Tie 
initial  height  of  HH  is  read  in  as  60  inches,  and  line  42  increases  HH 
by  1  inch  for  each  loop  of  a  "DO"  cycle.  For  each  value  of  0,  q,  HH, 
and  AN,  the  height  of  the  point  U  above  the  floor  is  called  ROOF(l), 
and  its  value  is  given  by 

R00F(I)  =  HH  cos  TCI)  +  [ AN ( K )  -  Q( I ) D  sin  T(I) 

where  HH  is  the  value  of  the  length  of  the  ordinate  NU. 

For  each  location  of  N  along  AB,  line  44  punches  out  the  first  value  of 
ROOF(l)  to  equal  or  exceed  76  inches.  Line  45  reduces  HH  to  74 
inches.  Next,  the  program  moves  N  4  inches  closer  to  B  and  repeats 
the  whole  process  until  N  is  at  the  end  B  of  the  underside.  Thus,  with 
the  108-inch  container  considered,  there  are  28  different  locations  for 
N.,  The  initial  value  of  AN,  namely,  AN(1),  which  is  read  in  for 
FORTRAN  Program  Number  Four,  is  -54  inches. 


In  the  line  of  data  following  the  "END"  or  Program  Number  Four,  the 
starting  station  for  the  forward  side  of  the  container  was  considered  to 
be  265  inches  from  the  bulkhead  at  the  forward  end  of  the  cargo  com¬ 
partment;  this  arrangement  corresponds  to  loading  plans  1  and  3.  A 
container  starting  this  close  to  L  has  the  greatest  period  of  time  to  tip 
and  thus  to  strike  the  ceiling.  Thus,  a  critical  case  is  being  used  in 
FORTRAN  Program  Number  Four;  nevertheless,  it  can  be  seen  that 
the  container  can  have  a  constant  height  of  76  inches  without  fouling 
the  ceiling.  The  same  applies  to  all  other  containers  starting  in  any  of 
the  six  loading  plans. 

NOTE:  As  soon  as  this  envelope  departs  for  any  reason  from  the 
simple  rectangular  outline  considered  in  this  report,  it  may  be 
necessary  to  adopt  a  more  sophisticated  expression  for  the  moment 
of  inertia  about  a  horizontal  axis  through  the  center  of  gravity.  A 
container  with  contents  of  nonuniform  density,  a  small  truck,  or  a 
gun  might  well  be  so  asymmetric  that  the  theory  set  forth  in  this 
study  might  require  modification. 


THREE-DEGREES- OF-FREEDOM  MOTION 


This  phase  occurs  after  the  container  leaves  the  floor.  Particular 
attention  is  paid  to  the  possibility  of  whether  the  inclined  ramp  at  the 
rear  of  the  floor  could  be  struck.  In  Figure  5,  axes  are  taken  through 


Di rection  of  Motion  of  Aircraft 


C 


the  end  L  of  the  floor --Lx 
positive  in  the  aft  direction 
and  Ly  positive  in  the  uown- 
ward  direction.  The  configu¬ 
ration  of  the  three-degree  sys¬ 
tem  is  defined  by  the  coordi¬ 
nates  x  and  y  of  the  center 
G  G 

of  gravity  and  by  0  ,  the  angle 
of  inclination  of  the  underside 
of  the  container  to  the  hori¬ 
zontal. 


The  potential  energy  of  the 
system  is  that  of  its  total 
mass  situated  at  the  center  of 
gravity  plus  the  work  done 

against  the  force  of  the  parachute  cord  from  some  standard  position  to 
its  present  position.  Thus, 


Figure  5.  Sketch  Illustrating  Three- 
Degrees -of- Freedom 
Motion. 
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rw  i  _  * 

V  =  constant  -  WyQ  -  P  cos  X(xQ  -  HG  sin  9  +  -5-  cos  0) 

AM 

+  P  sin  \(yc  t  HG  cos  0  ty  sin  6)  . 

The  kinetic  energy  of  the  system  is  that  of  its  total  mass  moving  with 
the  center  of  gravity  and  the  kinetic  energy  of  rotation  about  G.  Thus, 

2KEg  *  W(x§  +  y|)  +  IG02  . 

LAGRANGE  EQUATIONS  OF  MOTION  FOR  COORDINATES  x  AND  y 
AND  FOR  ANGLE  OF  INCLINATION 


For  the  coordinate  x, 

Wx„  -  Pg  cos  X  (6) 

G 

and  for  y, 

•• 

Wy„  =  g(W  -  P  sin  X) .  (7) 

For  the  angle  of  inclination,  0, 


AM 


I  *0  =  g(-p  HG  cos  0  cos  X  -  p  cos  X  sin  0 


AM 


+  P  HG  sin  X  sin  6  -  P  cos  0  sin  X)  .  (8) 


This  last  equation  yields 


IG6 


——  =  cos  0(-HG  cos  X  -  sin  X)  -  sin  9(-HG  sin  X  +  cos  X)  . 


AM 


gP 


2HQ 

Let  8  be  such  that  tan  3  =  then, 


10 


D  l— - — 

=-J hg2  +  (M)2  [-  cos  0  sin  (X  +  3)  -  sin  0  cos  (X  +  8)] 

W  1 


=  -~JhGz  +  (^r)  sin  (0  +  X  +  8). 
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The  square  of  the  radius  of  gyration  of  the  container  about  the  hori¬ 
zontal  axis  through  G  perpendicular  to  the  plane  of  the  paper  is 

AM2  +  4HG2 

- n -  . 


Hence, 


-9g  sin  (9  +  X  +  jO 
W4HG^  +  AM2 


(9) 


FORTRAN  Numerical  Solution  of  Equation  4>  =  -SS  sin 

Let  $*0+X+8.  In  the  FORTRAN  program,  ♦  is  recorded  as  F  and 
8  as  S.  Also  in  the  FORTRAN  program,  the  coefficient  of  sin  (0  +  X  +  8) 
in  the  right-hand  member  of  equation  (9)  is  called  -SS.  Therefore,  in 
the  FORTRAN  program,  the  exact  solution  of  $  =  -SS  sin  $  involves 
elliptic  integrals  of  the  first  kind,  which  are,  in  effect,  published  for 
discrete  values  of  SS.  To  find  4  for  intermediate  values  of  t  and  of  SS 
would  thus  involve  a  double  interpolation  by  t  and  SS  and  would  prove 
to  be  troublesome.  It  is  therefore  more  convenient  to  find  an  approxi¬ 
mate  numerical  solution  of  equation  (9)  in  the  same  way  as  has  been 
done  for  other  differential  equations. 

FORTRAN  Program  Number  Three  (see  appendix)  is  a  double  one:  it 
first  solves  equation  (9)  and  then  useB  the  values  of  $  and  0  to  obtain 
the  trajectory  of  corner  B  of  the  container  in  such  a  way  that  it  can  be 
readily  decided  whether  or  not  the  15-dcgree  ramp  is  struck.  This 
program  relates  to  a  108-inch-long  container  starting  from  the  rear¬ 
most  station  in  loading  plans  1  or  3,  and  this  is  a  critical  case  for 
striking  the  ramp.  In  this  program  the  angle  X,  or  YAL,  is  zero;  that 
is,  the  parachute  cord  has  been  assumed  to  be  horizontal.  (Nonzero 
angles  of  this  inclination  can  also  be  tested  without  difficulty.) 

Lines  1  through  17  of  this  program  provide  the  solution  of  equation  (9). 

HG  and  AM  are  read  in  as  usual.  T(l)  and  TD(1),  which  are  the  initial 
values  of  0  and  its  first  derivative  with  respect  to  time,  and  DD,  which 
is  the  arbitrary  interval  of  time  between  successive  configurations  of 
the  three-degree  system  and  which  is  taken  to  be  .005  second,  are 
also  read  in. 


There  is  a  FORTRAN  subroutine  for  ATAN  so  that  8,  whose  value  is 


ATAN 


can  be  found  with  no  difficulty.. 
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The  initial  values  of  T(l)  and  TD(1)  for  the  three-degree  motion  are  the 
final  values  for  the  corresponding  two-degree  motion.  In  the  case  of 
the  54-inch  container,  the  container  remains  horizontal  until  corner 
B  reaches  L  and  the  three-degree  motion  starts.  Thus,  the  small 
container  enjoys  no  two-degree  motion  and  the  initial  values  for  T(l) 
and  TD(1)  for  this  container  are  both  zero. 

As4>=0+X  +  8,  the  derivatives  of  $  and  of  9  with  respect  to  time  are 
equal.  This  is  the  explanation  of  line  10  of  FORTRAN  Program  Num¬ 
ber  Three.  Line  11  finds  FDD(l)  from  the  differential  equation  itself, 
while  line  12  differentiates  the  equation  and  uses  the  result  to  deter¬ 
mine  FDDD(l). 

Thus  F(l),  FD(1),  FDD(l),  and  FDDD(l)  are  quantities  that  can  be 
substituted  in  the  Taylor  series  for  F{2),  FD(2),  FDD(2),  FDDD(2), 
and  so  on,  successively,  for  any  value  of  index  "I"  in  the  FORTRAN 
program.  4>=0  +  X  +  Basa  function  of  time  having  been  found,  the 
corresponding  values  of  6  can  be  used  to  find  the  trajectory  of  corner 
B.  If  equation  (6)  is  integrated  twice  with  respect  to  time, 

*e  ■  <00  *  t<Ss>„  *  co.  X 

where  suffix  o  refers  to  the  initial  conditions.  Evidently, x^  initially 
equals 


AM  cos  0Q 


and  from  the  known  formula  for  distance  covered  at  constant  accelera¬ 
tion, 


<*G)0  =V3g(377  -  xst)  • 

Thus, 


xQ  =  cos  X  +  tVU59. 2(377  -  xst>  +  289. 2t2  (10) 

where  t  is  time  in  seconds  from  corner  B's  being  at  L. 

Similarly,  equation  (7)  yields 


=  'Vo  * 


t* 


1.5  sin  X) 


(11) 
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r-  n  -  -y  - ai'em. 


Based  on  these  coordinates  of  G,  those  of  B  are  as  follows: 

XB  *  XG  “  HG  S*n  ®  *  4|-  coa  ®» 
which,  in  view  of  equation  (10),  gives 

Xj  Sy  t  tVll59.2(377-xst)  +  28S.2t2  -  HG  sin  TH(I)  -  M  cos  TH(I) 
and 

yB  =  yQ  ♦  HG  cos  6  -  sin  0  *  -HG  +  t(yG)o 

♦  193. 2t2  ♦  HG  coa  TH(I)  -  iy  sin  TH(I). 

Here  TH(1)  means  the  value  of  0  at  configuration  index  "I".  The  symbol 
TH  is  used  simply  to  emphasize  that  the  three-degree  phase  is  being 
considered.  Either  of  the  two  following  methods  could  be  used  to  find 
the  initial  values  of  TH(1),  THD{1),  and  QD(1):  a  FORTRAN  subpro¬ 
gram  (which  would  interpolate  the  values  of  these  three  quantities  be¬ 
tween  the  two  successive  values  of  Q,  one  less  than  and  the  next  greater 
than  )  could  be  incorporated  in  Program  Number  One,  or  the  same 
interpolation  could  be  performed  by  hand.  However,  only  six  loading 
configurations  of  containers  in  the  cargo  compartment  are  being  con¬ 
sidered,  and  those  most  likely  to  lead  to  the  ramp's  being  struck  are 
loading  plans  1  and  3  with  the  108-inch  containers  and  loading  plans  2 
and  5  with  the  54-inch  containers.  Therefore,  there  are  only  two 
distinct  loading  plans  of  interest;  (y^^will  be  calculated  by  hand  and 
the  result  will  be  read  in  as  data. 

yg  -  q  sin  0  -  HG  cos  0. 

Hence, 

•  •  •  • 

yg  =  q  sin  0  +  q  cos  00  +  HG  sin  00, 

and 

•  AU 

(yG>o  8  QD(1)  sin  TH(1)  +~T  THD(1>  +  HG  TH(1)  THD(l)  (12) 

if  0  is  small. 

To  find  QD(1),  TH(1),  and  THD(l),  Program  Number  One  is  run  with 
data  corresponding  to  a  108-inch  container  starting  from  a  station  265 
inches  from  the  forward  bulkhead.  It  is  then  possible  to  find  the  two 
successive  entries  in  the  output  for  Q;  the  first,  Q(3),  is  less  than 
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108 


-  'i4  '«  r  ■v'.ir'.yi 


=  54  inches,  and  the  next,  Q(4),  is  greater  than  54  inches.  Follow¬ 
ing  are  the  values  of  Q(3)  and  Q(4),  together  with  the  corresponding 
values  of  T  and  the  first  differential  coefficients  of  the  two  coordinates, 
Q  and  T. 


Q(3)  =  52.092921 

QD(3)  =  360.  35773 
T{3)  =  . 18743423E-04 

TD{3)  =  .  37486 846E-U2 


Q(4)  =  55.724500 

QD(4)  =  365.95819 
T(4)  =  .  88831920E-04 

TD{4)  =  .  10269014E-01 


The  three -degree  motion  starts  when  Q  =  54  inches;  by  simple  inter¬ 
polation  for  tnis  value,  TH(1,  =  .  55549696E-04,  THD(l)  =  .  71727702E-02, 
and  QD(1)  =  363.  29876.  If  these  values  are  substituted  in  equation  (12), 

(y  )o  =  . 020x81  +  .387331  +  .000012  =  .407524 

When  this  value  is  now  substituted  in  equation  (11),  the  second  part  of 
FORTRAN  Program  Number  Three  can  be  constructed,  that  is,  lines 
19  through  27,  The  expressions  for  Xg  and  yg  are  built  up  piecemeal 
because  of  the  need  to  keep  the  lines  of  FORTRAN  down  to  72  charac¬ 
ters.  The  ratio  yg  :  Xg  is  calculated  as  shown  for  each  value  of  index 
"I";  it  does  not  equal  or  exceed  tan  15°  =  .  2679  for  any  Xg  less  than  88 
cos  15°  inches,  88  inches  being  the  length  of  the  ramp.  This  shows 
that  the  ramp  is  not  struck  in  loading  plans  1  or  3. 


Foi  loading  plans  2  and  5  with  the  54-inch  containers,  the  calculations 

are  more  simple.  Here,  AI\4  =  54  inches,  xgt  =  291  inches,  T(l)  =  0 

and  TD(1)  =  0,  Q(l)  =  27  inches,  and  QD(1)  =  315.739  inches  per 

second.  Thus,  (y  )  =  0,  and  the  first  term  in  line  24  of  FORTRAN 

Go 

Program  Number  Three  for  YBA(l)  drops  out.  Once  again  the  ratio 
yB:xB  never  equals  or  exceeds  tan  15°  for  Xg  less  than  88  cos  15 

inches,  so  it  can  be  concluded  that  the  ramp  is  not  struck. 


APPENDI X 


FORTRAN  PROGRAMS 

FORTRAN  PROGRAM  NUMBER  ONE 


C  SOLUTION  OF  TWO-DEGREES-OF-FREEDOM  EQUATIONS  i 

DIMENSION  T(21),TD(21),TDD(21),Q(21),QD(21),QDD(21)  2 

READ  3,T(l),TD(l),Q(l)tQD(l),AM,DE,YAL,HG  3 

C=386.4*((Q(l))*COS(T(l))-1.5*(Q(l)+AM/2.)*SIN(YAL))  4 

G=( ( ( AM**2 )+4 . *( HG**2 ) ) /12 . )+( HG**2 )  5 

B=HG  6 

A=G+(Q(1))**2  7 

F=386.4*(1.5)*(COS(YAL))t((TD(l))**2)*Q(l)  8 

E=l.  9 

D=HG  10 

TDD( 1 )=( C*E-F*B ) /( A*E-D*B )  11 

QDD(1)=(A*F-D*C)/(A*E-D*B)  12 

C  T  AND  Q  HAVE  NOW  BEEN  PREPARED  AT  POSITION  OF  SUBSCRIPT  (1)  13 

DO  51  1=1 „20  14 

Q( I+1)=Q( I )+DE*QD( I )+( ( ( DE**2 ) /2 . )*QDD( I ) )  15 

QD( I+1)=QD( I )+DE*QDD( I )  16 

T( I+1)=T( I )+DE*TD( I )+( ( DE**2 ) /2 . )*TDD( I )  17 

TD( I+1)=TD( I )+DE*TDD( I )  18 

CAA=0(I+1)*COS(T(I+1))+HG*SIN(T(I+1))  19 

CAB=-1.5*(Q(I+1)+AM/2.)*SIN(T(I+1)+YAL)  20 

CA=CAA+CAB  21 

CB=-2 .*( TD( 1+1 ) )*( QD( 1+1) )*Q( 1+1)  22 

CC=(386.4*CA)  +  CB  23 

BA=HG  24 

AA=G+(Q( 1+1) )**2  25 

FAA=386.4*(SIN(T( 1+1) )+l, 5*C0S( T( I+1)+YAL) )  2b 

FAB=( (TD( 1+1) )**2)*(Q( 1+1))  27 

FA=FAA  +  FAB  28 

EE=1.  29 

DD=HG  30 

TDD( 1+1) =( CC*EE-FA*BA) /( AA*EE-DD*BA )  31 

QDD(I+1)=(AA*FA-DD*CC)/(AA*EE-DD*BA)  32 

C  T  AND  Q  HAVE  NOW  BEEN  PREPARED  AT  ALL  POSITIONS  TO  SUBSCRIPT  21  33 

C  THERE  FOLLOW  TWO  FUNCTIONS  TO  TEST  THE  ACCURACY  OF  THE  DIFFL  EQUATS  34 
J=I+1  35 

FAMM=AA*TDD( J )+2 . *( Q( J ) )*( QD( J ) )*( TD( J ) )+( HG*( QDD( J ) )  )  36 

FAMN=386 . 4*(Q(J)*C0S(T( J ) )+HG*SIN(T( J ) ) )  37 

FAMP=579 , 6*(Q( J )+AM/2 , )*SIN(T( J )+YAL )  38 
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FAM*FAMM+FAMP-FAMN  39 

FAN A«HG*TDD( J ) - ( ( TD( J ) ) **2 ) *Q ( J ) +QDD( J )  40 

FANB»386,4*(SIN(T(J) )+1.5*C0S(T(J)+YAL) )  41 

FAN=FANA-FANB  42 

PUNCH  5,J,FAM,FAN  43 

PUNCH  4tJ,Q(I+l),QD(I+l),QDD(I+l)  44 

51  PUNCH  4 1 J ,T( I+l) ,TD( I+l) #TDD( I+l)  45 

STOP  46 

4  FORMAT  (I2,1X  E14.8,1X  E14.8.1X  E14.8)  47 

3  FORMAT  (F3.1,F3.1,F4.1,F5.1,F5.1,F3.2,F8.6fF4.l)  48 

5  FORMAT  (I3,1X  £14. 8, IX  E14.8)  49 

END  50 

0.00.037.7515.1216.0.010,03490720.0  51 


FORTRAN  PROGRAM  NUMBER  TWO 


C  INCLINATION  OF  SMALL  CONTAINER  DURING  FIRST  PHASE  1 

DIMENSION  X( 50 ) ,XD( 50 ) ,XDD( 50 ) ,T( 50 ) ,TD( 50 ) ,TDD( 50 )  2 

READ  3,AM,HGtDE  3 

READ  6,X(1),XD(1),T(1),TD(1)  4 

A=-HG  5 

Bs(AM**2  +  4.*HG**2)/3.  6 

C=-193.2*AM  7 

D=l.  8 

E=-HG  9 

Fs579.6  10 

XDD(1)=(C*E-F*B)/(A*E-D*B)  11 

TDD(l)s(A*F-D*C)/(A*E-D*B)  12 

DO  55  1*1,45  13 

X( I+l )=X( I )+DE*XD( I )+( ( DE**2 ) /2 . )*XDD( I )  14 

XD(I+1)*XD( I )+DE*XDD(I )  15 

T( I+1)*T( I )+DE*TD( I )+( ( DE**2 ) /2 . )*TDD( I )  16 

TD(I+1)*TD(I)  +  DE*TDD( I )  17 

AA=-((AM/2, ) )*SIN(T(I+1) )  +  HG*C0S(T(I+1))  18 

BB=(AM**2+4.*HG**2)/3.  19 

CCP=HG*SIN( T( I+l ))-( AM/2 . )*C0S(T(I+1) )-3,*(AM/2. )*SIN(T(I+1) )  20 

CC#386,4*CCP  21 

DD=1.  22 

EE=-( (AM/2, )*SIN(T(I+1) )+HG*CCS(T(I+l) ))  23 

FF=579, 6+((TD( I+l ) )**2)*(( AM/2 . )*C0S(T(I+1) )-HG*SIN(T(I+l)) )  24 

XDD( I+l )=( CC*EE-FF*BB ) /( AA*EE-DD*BB )  25 

TDD( I+l ) =( AA*FF-DD*CC ) /( AA*EE-DD*BB )  26 

PUNCH  4 , I ,X( I ) ,XD( I ) ,XDD( I )  27 

PUNCH  4,I,T(I),TD(I),TDD(I)  28 
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FORTRAN  PROGRAM  NUMBER  THREE 


Configurations  1  and  3  Trajectory  of  Edge  B  1 
Dimension  F( 51) ,FD( 51) ,FDD( 51) ,FDDD( 51) ,T(1) ,TD(1) ,TH(51)  2 
Dimension  TEE(  51) ,XBA( 51 ) ,XBB( 51 ) ,XB( 51 ) ,YBA( 51) ,YBB( 51) ,YB( 51)  3 
Dimension  RAT(51)  4 
Read  3,HG,AM,T( 1) ,TD( 1) ,DD  5 
Read  5,X  6 
SS=3477 . 6/( SQRT<  AM**2+4 . *HG**2 ) )  7 
S=ATAN(2.*HG/AM)  8 
F(1)=S  9 
FD(1)=TD(1)  10 
FDD(1)=-SS*SIN(T(1)+S)  11 
FDDD(1)=-SS*C0S(T(1)+S)*TD(1)  12 
DO  50  1=1,50  13 
F( I+1)=F( I )+DD*FD(  I )+( ( DD**2 )/2 . )*FDD( I )+( ( DD**3)/6 . )*FDDD( I )  14 
FD(I+1)=FD( I )+DD*FDD( I )+( ( DD**2 )/2 . )*FDDD( I )  15 
FDD(I+1)=-SS*SIN(T(I+1)+S)  16 
FDDD( I+l)=-SS*C0S( F( 1+1) )*FD( I+l)  17 
TH( I )=F( I )-S  18 
P=I  19 
TEE(I)=P*DD  20 
XBA(I)=TEE(I)*SQRT(1159.2*(377.-X))  r  289 .2*(TEE( I ) )**2  21 
XBB( I )=-HG*SIN( TH(  I ) )-( AM/2 . )*C0S(TH( I) )  22 
XB( I )=AM/2,+XBA(I )  +  XBB(I)  23 
YBA(I)=.4075*TEE(I)  +  193.2*(TEE(I))**2  24 
YBB( I )=HG*C0S(TH( I ) )  -  ( AM/2 . )*SIN(TH( I ) )  25 
YB(I)=-HG  +  YBA(I)  ♦  YBB(I)  26 
RAT( I )=  YB(I)/XB(I)  27 
PUNCH  4,1 ,XB(I ) ,YB(I ) ,RAT( I )  28 


50  CONTINUE 

STOP 

3  FORMAT  (F4.1,F5.1,E14.8,E14.8,F4.3) 

5  FORMAT  (F5.1) 

4  FORMAT  (13, IX  E14.8.1X  E14.8.1X  E14.8) 

END 

30.0108,0+.55549696E-04+.71727702E-02.005 

265.0 
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FORTRAN  PROGRAM  NUMBER  FOUR 


Lin* 


CALCULATION  OF  SAFE  ENVELOPE  1 

DIMENSION  T(21) fTD(21) tTDD(21) ,Q(21) ,QD{21) ,QDD(21) ,ROOF(21)  2 

DIMENSION  AN(30)  3 

READ  3,T(1) ,TD(1) ,X,AM,DE,YAL,HG,HH,AN(1)  4 

Q( 1)=( 1.5*( ( AM*SIN( YAL) )/2 ,+HG*COS(YAL) ) /( l.-1.5*SIN( YAL) ) )  5 

QD(1)*SQRT( 1159. 2*( 377 ,-X-( (AM/2 .-Q( 1) ) ) ) )  6 

C=386.4*((Q(1))*C0S(T(1))-1.5*(Q(1)+AM/2.)*SIN(YAL))  7 

G=( ( ( AM**2 )+4.*(HG**2 ) )/12 . )+( HG**2)  8 

B*HG  9 

A=G-f(Q(l))**2  10 

F=386.4*( 1.5*(C0S( YAL) ) )♦( (TD( 1) )**2)*Q( 1)  11 

E»l.  12 

D=HG  13 

TDD( 1 ) = ( C*E-F*B ) / ( A*E-D*B )  14 

QDD(1)=(A*F-D*C)/(A*E-D*B)  15 

39  DO  44  1=1,20  16 

Q( I+1)*Q( I )+DE*QD< I )+( ( ( DE**2 )/2 . )*QDD( I))  17 

QD(I+1)=QD(I)+DE*QDD(I)  18 

T( I+1)=T( I )+DE*TD( I )♦( ( ( DE**2)/2 . )*TDD( I))  19 

TD( 1+1 )=TD( I )+DE*TDD( I )  20 

CAA=Q(I+1)*C0S(T(I+1))+HG*(SIN(T(I+1)))  21 

CAB=-1.5*(Q(I+1)+AM/2.)*SIN(T(I+1)+YAL)  22 

CA=CAA+CAB  23 

CB=-2.*(TD(I+1))*(QD(I+1))*Q(I+1)  24 

CC=(386.4*CA)  +  CB  25 

BA=HG  26 

AA=G+(Q(I+1))**2  27 

FAA=386.4*(  SIN(  T(  1+1)  )+l .  5*C0S(  T(  I+D+YAL) )  28 

FAB=((TD(I+1) )**2)*(Q(I+1) )  29 

FA=FAA+FAB  30 

EE=] .  31 

DD=HG  32 

TDD( 1+1 )= ( CC*EE-FA*BA ) /( AA*EE-DD*BA )  33 

QDD( 1+1 )=( AA*FA-DD*CC ) /( AA*EE-DD*B A )  34 

44  CONTINUE  35 

DO  64  K=l,28  36 

P=K  37 

AN(K)=AN(1)  +  (P-1, )*4.  38 

38  DO  54  1=1,20  39 

ROOF( I )=HH*COS(T( I ) )+( AN(K )-Q( I ) )*SIN(T( I))  40 

IF  (ROOF(I)-76.)  54,55,55  41 

54  CONTINUE  42 

h'H=HH+l,  43 

GO  TO  38  j»4 

55  PUNCH  31 ,1 ,X ,ROOF( I ) ,HH , AN ( K )  45 

HH  =  74.  46 

64  CONTINUE  47 

STOP  48 

3  FORMAT  (F3,1,F3.1,F5.1,F5.1,F3.2,F8.6,F4.1,F4.1,F6.1)  49 

31  FORMAT  (13, IX  F5.0,1X  F9.4,1X  F5.0.1X  F5.0)  50 

END  51 


0.00.0265,0108.0.010.00000030.060.0-054.0 
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Wright -Patter son  AFB,  USA  Liaison  Det.  1 
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